Nebular Spectra of SN 1998bw Revisited: 
Detailed Study by One and Two Dimensional Models 
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ABSTRACT 

Refined one- and two-dimensional models for the nebular spectra of the hyper- energetic Type 
Ic supernova (SN) 1998bw, associated with the gamma-ray burst GRB980425, from 125 to 376 
days after B-band maximum are presented. One dimensional, spherically symmetric spectrum 
synthesis calculations show that reproducing features in the observed spectra, i.e., the sharply 
peaked [01] 6300A doublet and Mgl] 4570A emission, and the broad [Fell] blend around 5200A, 
requires the existence of a high-density 0-rich core expanding at low velocities (< 8,000 km 
s~^) and of Fe-rich material moving faster than the 0-rich material. Synthetic spectra at late 
phases from aspherical (bipolar) explosion models are also computed with a two-dimensional 
spectrum synthesis code. The above features are naturally explained by the aspherical model if 
the explosion is viewed from a direction close to the axis of symmetry 30°), since the aspherical 
model yields a high-density 0-rich region confined along the equatorial axis. By examining a large 
parameter space (in energy and mass), our best model gives following physical quantities: the 
kinetic energy £'51 = E^/W^^ ergs > 8 — 12 and the main-sequence mass of the progenitor star 
M^ms > 30 — 35Mq. The temporal spectral evolution of SN 1998bw also indicates mixing among 
Fe-, 0-, and C-rich regions, and highly clumpy structure. 



Subject headings: gamma rays: bursts - 
supernovae: individual (SN 1998bw) 



line: profiles - nuclear reactions, nucleosynthesis, abundances ■ 
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1. INTRODUCTION 

There has been accumulating evidence that 
a class of supernovae (SNe) is related to long- 
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duration Gamma-Ray Bursts (GRBs; see Piran 
1999 for a review) and possibly to their low energy 
analog, X-Ray Flashes (XRFs; Heise et al. 2001). 
The discovery of SN 1998bw in the error box of 
GRB 980425 (Galama et al. 1998) raised the is- 
sue, and now the association between GRBs and 
SNe is firmly confirmed by the emergence of su- 
pernova spectra in the optical afterglows of some 
GRBs, i.e., SN 2003dh/GRB 030329 (Hjorth et 
al. 2003; Kawabata et al. 2003; Matheson et al. 
2003; Stanek et al. 2003), SN 20031w/GRB 031203 
(Malesani et al. 2004; Thomsen et al. 2004, Maz- 
zali et al., in preparation), and SN 20021t/GRB 
021211 (Delia Valle et al. 2003). 

All these supernovae (except SN 20021t) seem to 
belong to a special class (see e.g., Matheson 2004 
for a review). Their spectra near maximum optical 
brightness are characterized by significant blend- 
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ing of very broad lines. A popular idea to explain 
this feature in the GRB-related supernovae is a 
hyper-energetic explosion (£^51 = ergs > 

5 - 10) of a very massive star (Af,„s > 20 - 25Mq). 
Such an energetic supernova is often called a "hy- 
pernova" (Iwamoto et al. 1998; Nomoto et al. 
2004) , with a somewhat different use of the origi- 
nal terminology suggested by Paczynski (1998) to 
describe the entire GRB/afterglow phenomenon. 

The nature of GRBs and hypernovae, and their 
mutual relation, are still the subject of debate. As- 
phericity seems to be the key to understand both 
GRBs and hypernovae. It is widely accepted that 
GRBs are produced by a relativistic jet viewed on 
or near its axis (Frail et al. 2001; Bloom, Frail, & 
Kulkarni 2003) . Their physical link to hypernovae 
then suggests that hypernovae also should show 
signatures of asphericity, according to popular 
progenitor scenarios: a black-hole plus an accre- 
tion disk (Woosley 1993; MacFadyen & Woosley 
1999; Brown et al. 2000), or a highly magne- 
tized neutron star (Nakamura 1998; Wheeler et 
al. 2000). 

The hypernova models for the GRB-related su- 
pernovae have been developed on the basis of mod- 
eling early phase observations up to 2 months 
after the explosion. For SN 1998bw, photometric 
and spectroscopic observations covering more than 
1 year are available (SoUerman et al. 2000; Patat 
et al. 2001). Since SN 1998bw was the first super- 
nova being identified as the counterpart of a GRB 
(GRB980425), and it has been extensively referred 
to as a prototypical hypernova, it is important to 
understand its nature. 

Studying the late phase light curve and spec- 
tra has provided additional hints on the nature 
of SN 1998bw. While the early phase observa- 
tions can be reproduced by a "spherical" hyper- 
nova model (Iwamoto et al. 1998; Woosley, East- 
man, & Schmidt 1999; Nakamura et al. 2001), 
the light curve after ^ 2 months declined less 
steeply than the spherical hypernova model pre- 
diction (McKenzie & Schaefer 1999). A similar 
problem exists for other hypernovae (e.g., Mazzali, 
Iwamoto, & Nomoto 2000; Yoshii et al. 2003) and 
some SNe Ilb/Ib/Ic (e.g., Clocchiatti & Wheeler 
1997). For hypernovae, Maeda et al. (2003a) at- 
tributed the failure of the spherical models to the 
possible existence of an inner low-velocity dense 
core, which may be formed as the consequence of 



an aspherical explosion (e.g., Maeda & Nomoto 
2003b). 

Late phase spectra provide an excellent tool to 
examine the geometry of the ejecta, as they can 
probe deeper than the early-phase spectra. The 
nebular spectra of SN 1998bw showed peculiar fea- 
tures which were not explained by a simple picture 
(Patat et al. 2001). First, the [01] 6300A & 6363A 
doublet shows a sharply peaked emission profile. 
The profile is different from both the parabolic 
profile expected from a homogeneous spherical dis- 
tribution of emitting material and the flat-topped 
boxy profile expected from material distributed in 
a shell. Sollerman et al. (2000) found that both 
the spherical hypernova model C0138 of Iwamoto 
et al. (1998) (CO138E30; the explosion with 
E51 = 30 of a I3.8M0 CO star evolved from a 
zero-age main sequence star of Mms ~ 4OM0), and 
C06C of Woosley et al. (1999) {E^i = 22, a 6.55 
Mq go star, and M^s ~ 25M0) yielded too broad 
[01] 6300A emission. The feature near 52OOA, 
interpreted as a blend of [Fell] in Mazzali et al. 
(2001) and Maeda et al. (2002) also shows a pe- 
culiarity: it is broader than the [01] 6300A. Maz- 
zali et al. (2001) showed that each line contribut- 
ing to the [Fell] blend should be broader than 
the [01] 6300 A, i.e.. velocities of once-ionized iron 
ions should be higher than those of neutral oxy- 
gen along the line of sight. Subsequently, Maeda 
et al. (2002) computed the theoretical line profiles 
expected from bipolar aspherical models. They 
showed that the peculiarity of the [01] 6300A and 
[Fell] 52OOA emissions in SN 1998bw can be ex- 
plained by the aspherical explosion models, if the 
axis of symmetry is roughly aligned with the line 
of sight. 

Although these previous studies revealed the 
peculiarities and suggested possible solutions, 
there were still some limitations. As for one di- 
mensional models, previous studies (Sollerman et 
al. 2000; Mazzah et al. 2001) did not model the 
full set of late-phase spectra at various epochs on 
the basis of a hydrodynamic model. Concerning 
two dimensional models, the analysis of Maeda 
et al. (2002) had the following limitations. (1) 
First, they applied a simplified time- independent 
computation for the line profiles. (2) They did not 
model the flux. A more detailed analysis based 
on nebular line emission physics is necessary to 
derive the nature, e.g., the explosion energy, of 
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the explosion. (3) Also, it should be examined 
if their models, focusing on the [Fell] 5200A and 
the [01] 6300A, are consistent with other emis- 
sion Hncs, e.g., Mgl] 4570 A and [Call] 7300 A. 
(4) Their model grid (in energy and mass) was 
rather coarse. To constrain physical values such 
as the energy, models with various energies and 
masses should be examined comprehensively. (5) 
Finally, they only examined one nebular spectrum 
(+216 days after the maximum optical luminos- 
ity). It is important to examine if "one" model 
can consistently reproduce the temporal evolution 
of spectra. 

The purpose of this paper is to perform a 
more detailed analysis of the nebular spectra of 
SN 1998bw than previous works. In addition 
to one dimensional spectrum synthesis, we use a 
two dimensional nebular spectrum synthesis code, 
which includes two-dimensional 7-ray transfer and 
a computation of ionization/thermal structures for 
two-dimensional density and abundance distribu- 
tions. First, the nebular spectra of SN 1998bw are 
examined by means of one dimensional spectrum 
synthesis in §2, to obtain conditions necessary to 
explain the spectra of SN 1998bw. The subsequent 
parts are devoted to the two dimensional mod- 
els. In §3 we describe models and method, and 
present a comparison between the synthetic spec- 
tra and the observations. Summary and discussion 
are presented in §4. 

2. ONE DIMENSIONAL MODELS 

2.1. Method 

We use a non-LTE code to compute nebular 
spectra. The code is an extension of the one-zone 
nebular code of Mazzali et al. (2001), allowing 
radially stratified density and abundance distri- 
butions as an input model. First, a Monte-Carlo 
Method (e.g., Cappellaro et al. 1997; Chugai 2000; 
Maeda et al. 2003a) is used to solve the trans- 
port of 7-rays released by the decay chain ^^Ni — » 
^^Co^ ^^Fe. Positrons released by the same decay 
chain are assumed to be trapped in situ (Axelrod 
1980). The light curve of SN 1998bw suggests that 
at a few hundred days 7-rays are still the dominant 
heating source (Nakamura et al. 2001; Maeda et 
al. 2003a), making the heating of the ejecta little 
sensitive to the detail of positron transport (e.g., 
Milne, The, & Leising 2001). 



Next, ionization and NLTE thermal balance in 
each shell are solved on the basis of the prescrip- 
tion given by Ruiz-Lapuente & Lucy (1992). For 
the ionization source, we consider only impact ion- 
izations by high-energy particles produced by the 
deposition of 7-rays (and e+). The ionization bal- 
ance is then solved by equating the impact ion- 
ization rate to the recombination rate for each 
ion. An assumption here is that photoionization is 
negligible. Although it is still a matter of debate 
(see e.g., SoUerman et al. 2004), it is argued that 
large multiple resonance scattering at high ener- 
gies effectively reduces the UV radiation, therefore 
reduces the photoionizations (Kozma & Fransson 
1998b; Kozma et al. 2005). Level populations are 
obtained by solving rate equations in steady state 
for each level along with the thermal balance (i.e., 
equating of the non-thermal heating rate and the 
line cooling rate). We consider mainly forbidden 
lines as a source of cooling (Ruiz-Lapuente & Lucy 
1992), but strong allowed transitions, e.g., Call 
IR, are also included. Under nebular conditions, 
radiative losses are entirely due to coUisionally ex- 
cited lines. We do not take radiation transfer ef- 
fects into account, since the epochs examined in 
the present work are later than 100 days so that 
radiation transfer effects are usually small at least 
at optical wavelengths (sec §4.5.1 for more detail). 
Observationally, after day 100 ~ 200, neither the 
feature at 4570A nor that at 6300A evolved signif- 
icantly (Patat et al. 2001), indicating that these 
features are totally dominated by forbidden lines 
and that radiation transfer effects are negligible. 
In particular, we do not include line scattering and 
fluorescence, an assumption that is usually made 
in nebular analysis (e.g., Kozma et al. 2005). Flu- 
orescence may have effects in some lines, i.e., Ca 
II H, K and IR-triplet, and Nal (Kozma & Frans- 
son 1998b), for which we do not present detailed 
modeling in the present work. Finally, a synthetic 
spectrum is obtained by integrating the emission 
from each shell. 

Throughout the paper, synthetic spectra are 
compared with the spectra of SN 1998bw taken 
at 125, 200, 337, and 376 days after B maximum 
(Patat et al. 2001). The first spectrum (+125^^) 
may not fully be nebular and hence can not be 
fully interpreted using the present computational 
code (see §4). Despite this, we include this spec- 
trum in the comparison. The spectrum is at least 
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partly nebular and therefore can still be used to 
reject any model producing a spectrum very dif- 
ferent from it. 

2.2. SN 1998bw: One Dimensional Ap- 
proach 

With the time-dependent computational code, 

we first examine the synthetic spectra for the 
original spherical hypernova model CO138E30 
(Iwamotoetal. 1998; Nakamura et al. 2001). The 
mass cut, i.e., the boundary between the collaps- 
ing core and the ejecta, chosen in the present work 
is deeper than the original, yielding the mass of 
56]\jj ^ O.SMq in order to roughly fit luminosities 
at the late phases (The original model with ^^Ni 
~ 0.4 — O.6M0 is fainter than the late-phase ob- 
servations; Nakamura ct al. 2001). Figure 1 shows 
the temporal evolution of the synthetic spectra as 
compared with the observed ones. We next try to 
obtain better fits by changing the density struc- 
ture below V < 10, 000 km s~^ and the abundance 
distribution throughout the ejecta (described be- 
low in more detail). The temporal sequence of 
the spectra of this "modified" CO138E30 model 
is shown in Figure 2. The density structure and 
abundance distribution are shown in Figures 3 and 
4, respectively. 

The "original" model has difficulties in repro- 
ducing both (a) line profiles in each single spec- 
trum and (b) temporal evolution. As for the line 
profiles (a), elements in outer layers such as O 
and Mg produce synthetic lines with very wide, 
flat-topped profiles, while the observed [01] 6300 A 
and Mgl] 4570A are narrow and sharply peaked. 
The flat-topped profile is the consequence of the 
distribution of emitting material with a central 
hole (See Figures 3 and 4a: O is distributed at 
> 12,000 km s~^). Contrary to the observa- 
tions, the synthetic [Fell] 5200A blend is narrower 
than the synthetic [01] 6300A. This is a typical 
character of any spherically symmetric evolution- 
collapse-explosion model, as was pointed out by 
Mazzali et al. (2001). The problem in the tempo- 
ral evolution (b) is that the "original" CO138E30 
model yields very rapidly fading lines [01] 6300A 
and Mgl] 4570A between the epochs +125'^ and 
-^200^^ (Fig. 1). 

The modified model (Fig. 2) is constructed so 
as to overcome these problems. The line profiles 



suggest that Fc (mainly synthesized as ^^Ni at the 
explosion) is on average distributed at higher ve- 
locities than hydrostatic burning products, e.g., O 
and Mg. The sharply peaked [01] and Mgl] lines 
can be accounted for if these lines are predomi- 
nantly emitted at low velocity. In the one dimen- 
sional representation, a peculiar abundance distri- 
bution is therefore required: A high density 0-rich 
core is added to the original CO138E30 model. It 
also helps to explain the more rapid decline of the 
[Fell] than the [01] 6300A and the Mgl] 4570A 
emission (Patat et al. 2001), since the contribu- 
tion of higher velocity material becomes smaller 
as time goes by. Because deposition rate is larger 
in the modified model than in the original one, 
the mass of ^^Ni is smaller in the former, M(^^Ni) 
~ 0.6Mq. 

Another feature of the modified model is some 
mixing between the Fe-rich layer (originally at ve- 
locities at ~ 8, 000-12, OOO km s'^) and the 0-rich 
core (added at < 8,000 km s~^). The observed 
[Fell] 52OOA and [Call] 7300A fines show mildly 
peaked profiles (though less sharply peaked than 
the [01] and Mgl]), suggesting that a small frac- 
tion of the Fe-rich materials are mixed downward 
into the low velocity 0-rich core. The suggestion 
that (at least in a spherically symmetric represen- 
tation) mixing is necessary was already made by 
SoUerman et al. (2000). Our model represents 
an even more unusual structure (from the point of 
view of ID hydrodynamics), since the 0-rich core 
is located at low velocities together with just a 
small fraction of Fe. This helps to keep the ioniza- 
tion of Fe as high as what is inferred from the ob- 
served spectra, although the usual mixing (SoUer- 
man et al. 2000) with an Fe-dominated central 
region will lead to too low an ionization. 

The mixing is also favored in view of the tem- 
poral evolution of each line intensity. Without any 
mixing, energy deposited by 7-rays and positrons 
into the 0-rich region is all reprocessed into oxy- 
gen lines, predominantly into the [01] 6300, 6363A 
doublet. Therefore, the temporal evolution of [01] 
6300A emission roughly traces the deposition rate, 
yielding too rapid fading of the line from -1-125*^ 
to -1-200'* (Fig- !)• If other elements such as iron 
are mixed in, then the situation is different. In- 
tensity of each line depends on contribution from 
others, which also depends on thermal condition 
within the nebula. In the materials composed with 
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both iron and oxygen, the intensity ratio [Fell] 
5200 A/ [01] 6300 A decreases rapidly with decreas- 
ing density, yielding larger relative intensity of [01] 
6300A in more advanced epochs. Given that the 
sum of the intensities of these lines roughly follows 
the deposition rate, increasing relative importance 
of [01] 6300A with time yields the [01] loss rapidly 
fading than in the absence of any mixing. 

The modified model (Figs. 3 and 4b) is con- 
structed by mixing Fe from the Fe-rich layer down- 
ward into the 0-rich core and compensate for it by 
mixing oxygen from the 0-rich core into the Fe- 
rich layer. The same is also done for Ca. In addi- 
tion, intensity of the [01] 6300A at the first epoch 
is further decreased by introducing clumping in 
the 0-rich region. Regarding clumping, first in- 
troduced for SN 1998bw by Mazzah et al. (2001), 
the following facts give hints: (1) Without clump- 
ing [Felll] is too strong at ^ 47GGA, which is not 
seen in the observations. (2) The intensity ratio of 
[Call] 7300A to Call IR is too large in the original 
CO138E30 model. Therefore, we have introduced 
a filling factor of 0.1 throughout the ejecta The 
filling factor is kept unchanged at all epochs. 

Comparing Figures 1 and 2, it is seen that the 
fit to the observed spectral sequence is very much 
improved. There are still some difficulty in repro- 
ducing the spectrum at -t-125'^ (while it is much 
better than the original model). In addition to 
the strong [01] 6300A emission in the synthetic 
spectrum, Mgl] 4570A was also too strong at the 
first epoch (In Fig. 2, the mass fraction of Mg is 
reduced by a factor of 10 at the first epoch rela- 
tive to the subsequent epochs). However, it may 
simply be that the density in the first epoch is 
too high to apply the nebular spectral computa- 
tion (see §4). The present model is in a sense 
ad hoc, introducing peculiar elements distribution 
and density structure (i.e., 0-rich core). We will 
show in §3 that this apparent peculiarity in one- 
dimensional representation can naturally be inter- 
preted in the context of a two-dimensional model. 

3. TWO DIMENSIONAL MODELS 

3.1. Method and Models 

We have developed a two-dimensional code, 
which is applicable to any two-dimensional distri- 
bution of density and abundances (including the 
distribution of the heating source ^^Ni). In the 



current version, axisymmetry along the z-axis (po- 
lar axis) and reflection symmetry on the equatorial 
plane are assumed. The included physics is the 
same with the one dimensional code. The differ- 
ence between the one- and two-dimensional ver- 
sions is in the treatment of the 7-ray deposition 
and the computation of the whole spectrum. The 
7-ray deposition is solved as a two dimensional ra- 
diation transport problem using the Monte-Carlo 
Method. Ionization and NLTE thermal balance 
are solved in each mesh zone independently, since 
these processes take place locally as long as the op- 
tical depth is negligibly small. A synthetic spec- 
trum is then computed for different orientations 
(divided uniformly into 10 degree angular zones), 
taking into account different Doppler shifts for dif- 
ferent orientations. 

We use as input the models presented in Maeda 
et al. (2002). These include spherical and aspheri- 
cal explosions of a 16Mq He star, the core of a star 
of Mms = 4OM0. Asphericity was generated as- 
suming angle-dependent energy deposition, pref- 
erentially concentrated toward the polar direction, 
at the center of the collapsing core (See Maeda et 
al. (2002) for details). 

The models we explore in this paper are listed 
in Table 1. The distributions of density and of 
a few selected elements in the homologous expan- 
sion phase are presented in Figure 5. We examine 
models with various degree of asphericity defined 
by the parameter BP. The value of BP is the axial 
ratio of the initial aspherical energy injection. See 
Maeda et al. (2002) for details (their parameter 
a//3 is the same as BP). The models correspond 
to those of Maeda et al. (2002) as follows: BP16 
(A), BPS (C), BP2 (E), and a spherical model BPl 
(F). In the present work, we add an additional 
model (BP4) which has an asphca'icity intermedi- 
ate between models BPS and BP2, by repeating 
the calculations of Maeda et al. (2002) but with a 
different asphericity parameter. 

As seen in Table 1, we generate a series of mod- 
els with various kinetic energy by artificially multi- 
plying the velocities at any points in the ejecta by 
a factor /y. The factor fv is taken as a parameter, 
ranging from 0.7 to 1.9. For each model specified 
by the parameters BP and fy, the masses of oxy- 
gen (Mo) and ^^Ni in the ejecta are changed by 
hand in order to obtain a good fit to the intensities 
of the [01] 6300A, 6363A and the [Fell] 5200A (see 
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§3.2 for details). The kinetic energy in the cjecta 
is scaled as Ek oc Mejfv^, where Mej is the ejecta 
mass, mainly oxygen, obtained by spectral fitting. 
Because of the procedure, the oxygen mass and 
therefore the kinetic energy are different from the 
original values (Mq ^ 8Mq and E51 ~ 10). 

Computing nebular spectra for these models, 
we obtain synthetic spectra resulting from various 
energies and cjecta masses (and thus progenitor 
masses). Summarizing, a spectrtun of a partic- 
ular model is specified by three parameters: the 
asphericity BP, the velocity scale fy (which gives 
the energy in the ejecta), and the orientation to 
the observer 6 (defined as the angle between the 
polar and the observer direction). 

3.2. Overall Spectra and Model Construc- 
tion 

We compare our synthetic spectra for the two 
dimensional, aspherical supernovac with the spec- 
tra of SN 1998bw (Patat et al. 2001) obtained at 
125, 200, 337, and 376 days after B maximum. We 
especially focus on the +337"^ spectrum, since the 
epoch is sufficiently late that the contribution of 
allowed Fe transitions (which are not included in 
the model) on the [Fell] blend near 5200A is likely 
negligible (see §4.5.1). 

Figure 6 shows examples of the synthetic spec- 
tra with a viewing angle 9 = 30°. The synthetic 
spectra are very much improved for the aspherical 
model, especially for the one with BP=8 and fv 
=1.15, compared to the spherical hydrodynamic 
explosion models (either the original CO138E30 
in Fig. 1 or models BPl in Fig. 6). The sharply 
peaked Mgl] 4570A, [01] 6300, 6363A, and Nal 
59OOA, as well as the broad [Fell], are naturally 
explained. Indeed, the fit by the aspherical model 
BPS (/y=1.15) is as good as the one by the one 
dimensional "modified" CO138E30 model. The 
masses of oxygen and ^^Ni are determined so as 
to obtain the correct line strengths (at -|-337days) 
in the two-dimensional models. This is done by re- 
ducing/increasing the density uniformly through- 
out either the ^^Ni or O rich region (see below). 
Therefore, each line profile is not sensitively af- 
fected by this procedure but results from the ex- 
plosion model itself. 

Figure 6 also shows how different kinetic en- 
ergy and asphericity affect the synthetic spectrum. 



For example, we can definitely rule out the model 
BP2 with /v=1.6 for SN 1998bw, since the [01] 
6300A line is too broad and its shape is less sharply 
peaked than in the observation. We discuss the 
profiles of each line in more detail in §3.3 and in 
the Appendix. The selection of models acceptable 
for SN 1998bw is discussed in §3.4. 

Because the intensities of the [Fell] 5200A and 
the [01] 6300A are well defined in the observation, 
they can be used to constrain the masses of iron 
(mostly the product of ^^Ni decay) and oxygen. 
These are model-dependent as the combination of 
mass and energy gives the average density, which 
in turn gives the deposition rate. For example, 
we find that the original hydrodynamic models 
(where Mq ~ 8Mq and S51 ~ 10) yield too strong 
[01] emission relative to [Fell] 5200A. Therefore 
Mo should be reduced to fit the line ratio cor- 
rectly (e.g.. Mo ~ 4.7Mq in model BPl, /y = l). 
In the current study, we have uniformly reduced 
densities in oxygen-dominated regions, and there- 
fore Mq, until we obtained the correct intensity 
ratio [01] 6300A/[FeII] 5200A at day -h337. The 
mass of ^^Ni is also constrained by the total lu- 
minosity in the observed spectrum. Therefore, we 
also changed the mass of ^^Ni (densities in ^^Ni 
dominated regions) until we obtained the correct 
luminosity at -^337'^. The masses M^^^Ni) and 
Mo consistent with the observed flux at -1-337'^ 
are thus derived for all the models as listed in Ta- 
ble 1. Note that a model with larger fv has lower 
density and smaller deposition rate, and therefore 
needs larger Mq and M(^^Ni) to fit the observed 
intensities. This is the reason the energy Ek in 
the larger fv models is larger than expected from 
the simple scaling Ek oc /y^, because the mass of 
the ejecta needed is also larger for larger fv- 

3.3. Emission Lines 

Dependencies of line profiles of the [Fell] blend 
52OOA, the [01] 6300, 6363A doublet, the Mgl] 
4570 A, and the [Call] 7291, 7324A on the three 
parameters, fv, 0, and BP are described in Ap- 
pendix in detail. Here we give a very brief sum- 
mary. 

• fv'- Irrespective of BP and 6, larger fv 
yields broader profiles for every line. 

• 6: The profile depends on whether the emit- 
ting element is explosive burning product 
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(Fc and Ca) or hydrostatic product (O and 
Mg). For the former, the hne is broader for 
smaller 6. For the latter, the dependence is 
in the opposite sense. 

• BP: For explosive burning products, larger 
BP yields a broader [narrower] line if 9 is 
small [large]. For hydrostatic products, the 
dependence is in the opposite sense. 

3.4. SN 1998bw: Two Dimensional Ap- 
proach 

In this section, wc discuss which models are ac- 
ceptable for SN 1998bw, according to each line 
profile (§3.3 and Appendix). On the basis of the 
selected models, we further seek for conditions 
that may yield better fits to the observations in a 
manner similar to what was done for the ID model 
(by introducing mixing and clumping; §2.2). 

Figure 7 illustrates how the model parame- 
ters are constrained by the observations. A fit 
is judged by the following rules. For the [Fell] 
5200A. it is checked if the width of the blend is 
consistent with the observed value. The region of 
acceptance is larger for ^ = 0° than 6 = 30°, since 
the former results in a broader profile. For the 
[01] 6300A, if a model fits well the line wings and 
does not produce a double-peaked emission line, it 
is regarded to be acceptable. The fit to the Mgl] 
4570A is not included, because the Mgl] gives the 
same constraints as the [01] 6300A lines, but with 
less accuracy. 

The fit to the [Call] 7300A line is uncertain. In 
a strict sense, no present model, either spherical or 
asphcrical, gives acceptable fits, since the observed 
line center is shifted blueward relative to the mod- 
els. See Appendix for details. Here, we tentatively 
regard a fit as acceptable if the emission at the 
blue wing is consistent with the observation. By 
doing this, models producing a doubled-peaked or 
a too broad flat-topped profile are excluded. There 
are several ways to shift the line center or to de- 
press the emission in the red, as we will discuss 
in §4. The region of acceptance for [Call] 7300 A 
moves to more aspherical and more energetic mod- 
els (yielding a broader Ca line) for larger viewing 
angle 6 (yielding a narrower Ca line), so that these 
two effects compensate one another. 

Figure 7 shows that both the [Fell] 5200A and 
the [01] 6300A are explained by highly aspherical 



(BP > 8), energetic models [fv > 1) viewed in 
a direction close to the pole {9 < 30°). Viewing 
angles 6 » 30° do not produce an acceptable fit 
to [01] 6300A. Also, the width of the [01] Hne 
sets an upper limit for fv (and therefore on the 
energy), since models with fv > 1-6 produce too 
broad [01] emission. If we include [Call] 7300A 
in the fit despite the large uncertainty, 9 ~ 30° 
is preferred, since highly aspherical and energetic 
models produce double peaked or too broad [Call] 
emission if viewed at ^ ~ 0° or even at ^ ~ 15° 
as mentioned in Maeda et al. (2002) (who did not 
include the [Call] in the model) . If only the [Fell] , 
[01] (and Mgl]) are considered in the fit, the fit by 
models with ^ = 15° is as nice as ^ = 30°. The 
[Call] 7300A feature is complicated and therefore 
needs further study to understand its nature (§4). 

Next, we examine the temporal spectral evo- 
lution of one of the acceptable models (BPS, 
/y=1.15, and 6 = 30°). This is shown in Fig- 
ure 8a. Additionally, Figure 8b shows model BP8 
with /y=1.6 and 9 = 30°. The latter model is not 
regarded as "acceptable" because the [01] 6300A 
line is too broad if it is normalized at the emission 
peak (see Appendix). The fit to the wings of the 
[01] emission, however, is not bad, therefore the 
fit is marginal. 

For SN 1998bw, this is the first attempt to 
check the consistency of a theoretical model by 
time-dependent computations. Figure 8 shows 
that both models have difficulty in reproducing 
the slowly fading [01] 6300A and Mgl] 4570A from 
-1-125"^ to -|-200^, a similar problem also encoun- 
tered for the one dimensional model CO138E30. 

As done in §2, we have attempted to obtain 
better fits, by introducing mixing and clumping. 
Guided by the ID model, we have tried to simulate 
these effects. First, mixing between the Fe-rich re- 
gion and the 0-rich region is taken into account by 
changing a fraction of O in the 0-rich region into 
Fe and Ca. Second, wc have also increased the 
mass of C in the 0-rich region since the synthetic 
feature near 8500A shows a deficiency of [CI] (Fig. 
8). Obviously liydrodynamic explosion models, ei- 
ther one or two dimensional, yield a mass fraction 
of carbon in the 0-rich region insufficiently small 
to fit the [CI] line strength (sec also Fig. 1). Fi- 
nally, we have introduced filling factors to simu- 
late clumping both for the Fe-rich and the 0-rich 
regions. 
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These parameters for mixing and clumping are 
set so as to obtain a fit as nice as possible. Figure 
9 shows the spectra of the "modified" 2D models. 
The parameters for the modification are as follows. 
For the BPS, /y=1.15 model, 3%, 1%, and 10% 
of O in the 0-rich region are changed to Fe, Ca, 
and C, respectively. Then filling factors of 0.1 and 
0.2 are introduced in the Fe- and 0-rich regions, 
respectively. For the BPS, /y=1.6 model, 1.5%, 
1%, and 10% of O in the 0-rich region arc changed 
to Fe, Ca, and C, respectively. Then filling factors 
of 0.05 and 0.1 is introduced in the Fe- and 0-rich 
regions, respectively. 

Comparing Figures 8 and 9 shows that the fit 
is improved very much by introducing mixing and 
clumping. The modified 2D models match the ob- 
served spectra as nicely as the modified ID model 
(see Fig. 9). The main deference between the 
modified ID model and the 2D model is that the 
modified ID model is very different from the orig- 
inal CO13SE30 model. The modification yielded 
line profiles different from the original model. The 
distributions of density and abundances in this 
model are inconsistent with the results of "ID" 
explosion hydrodynamics. In the 2D model, on 
the other hand, the distribution of density and 
abundances are based on the aspherical explosion 
model. The changes due to the modification of 
the abundance distribution and the introduction 
of mixing and clumping, which may naturally oc- 
cur, are small. This does not sensitively affect the 
line profiles. 

4. Summary and Discussion 
4.1. Summary 

We have revisited the nebular spectra of SN 

199Sbw, using one- and two-dimensional nebular 
spectrum synthesis codes. Compared with pre- 
vious works (ID models by SoUerman et al. 2000 
and Mazzali et al. 2001, and a 2D model by Maeda 
et al. 2002), we have extended the analysis in the 
following points. In the ID model, we have tried 
to construct a model which explains the temporal 
evolution of the spectra of SN 1998bw from +125'* 
to -1-376'^. The distribution of the density and ele- 
ments of the model (the modified CO138E30) is far 
from what is expected from ID spherically sym- 
metric explosion models. We have further com- 
puted 2D synthetic spectra of the aspherical explo- 



sion models of Maeda et al. (2002). These models 
show that an aspherical explosion scenario natu- 
rally explains the late-phase spectra of SN 1998bw. 
In the 2D model, we have extended the analysis 
of Maeda et al. (2002), including more realistic 
computations and more detailed comparison with 
the observations. 

As for the one- dimensional model, our main re- 
sults are as follows. (1) The original CO13SE30 
model for SN 1998bw (obtained by fitting the early 
phase observations up to ~ 2 months) is not con- 
sistent with the late-phase spectra. (2) A high- 
density 0-rich core below ~ 8, 000 km s^^ is nec- 
essary to explain the sharply-peaked [01] 6300A, 
6363A doublet profile. (3) At the same time, 
a large amount of ^^Ni (which decays into Fe) 
should be located at u > 8,000 km s~* to repro- 
duce the large ratio between the widths of [Fell] 
5200A and [01] 6300A. The above three features 
are not accounted for by any spherically symmet- 
ric evolution-collapse-explosion calculations. (4) 
Heavy elements such as Fe and Ca are likely mixed 
into the 0-rich region, and vice versa. Namely, O, 
Fe, and Ca do not occupy completely separated 
layers. (5) Clumping is important both in 0-rich 
and Fe-rich regions. 

In the context of the two-dimensional models, 
our results can be summarized as follows. (1) The 
[Fell] 5200A and the [01] 6300A lines are naturally 
reproduced by (some) aspherical models, confirm- 
ing the suggestion of Maeda et al. (2002). (2) 
The profile of the Mgl] 4571 A gives additional sup- 
port to the aspherical model for SN 1998bw. The 
analysis of the [Call] profile demonstrates, how- 
ever, that the present models are still not perfect. 

(3) Examining a wide parameter range, we have 
shown that nebular spectra of SN 1998bw indeed 
require a hyper-encrgctic explosion (sec also §4.2). 

(4) Similar to the ID model, mixing and clumping 
are at least partly responsible for the slow spec- 
trum evolution from -hl25^ to -1-200''. 

4.2. The Nature of SN 1998bw 

According to the present model, we can infer 
the nature of SN 1998bw. It was an highly aspher- 
ical explosion (BP ~ 8) viewed at 0° < 6* < 30°. 
The viewing angle is consistent with the off-axis 
model (16° < ^ < 36°) explaing properties of 
prompt gamma ray emission of GRB 980425 as- 
sociated with SN 1998bw (Yamazali, Yonetoku, & 
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Nakamura 2003) . The kinetic energy in the explo- 
sion models giving the best fits was £^51 ~ 8 — 12, 
mass of 5^Ni M{^'^Wi) ~ 0.4 - 0.5Mq), and mass 
of oxygen in the cjccta Mq ^ 5 — 6Mq. The mass 
of oxygen would correspond to a CO core of ~ 
9 - IOM0 and to a He core of ~ 11 - 12Mq, which 
is evolved from a star with My^s ^ 30 — 35M0. 
The argument assumes the central remnant is a 
AMq black hole, which could be larger if a larger 
fraction of the 0-dominated layer is accreted onto 
the black hole (e.g., MacFadyen & Woosley 1999; 
Maeda & Nomoto 2003b). 

Here we would like to discuss the uniqueness of 
the model. A fit to each line is used to restrict 
a possible range in the parameter space (BP, fy, 
and 0). Because the line shapes of different ele- 
ments depend on the parameters differently (see 
§3.3 and §3.4), the combination of fits to sev- 
eral lines is useful to narrow down the parameter 
space (Fig. 7). Once these parameters are given, 
then the masses of ^^Ni and ejecta mass are de- 
rived based on the line intensities, and therefore 
M(^^Ni), Mejectai and Ek are determined rather 
uniquely. The above model is derived including a 
qualitative fit to [CaH] 7300A, which is not very 
certain in the present study (see §3.4 and §4.5.2). 
If we omit this line, and use only the [01] 6300A, 
Mgl] 4570A, and [Fell] 5200A, then the degener- 
acy is not completely resolved (Fig. 7). For ex- 
ample, assuming 9 = 30° the models with BP > 8 
and 1 < /y < 1.45 (6 < £^51 < 16) are accept- 
able. For smaller 9, the range is even larger, while 
values ^ ~ 0° are disfavored because they do not 
even give a qualitatively acceptable fit to [Call] 
73OOA. The models can be further constrained by 
fitting early-phase observations, i.e., light curve 
and spectra (K. Macda, in preparation). 

We regard the above mass of oxygen and es- 
pecially the kinetic energy in the two-dimensional 
model as lower limits of those in the real ejecta for 
the following reason. Although we have computed 
a spectrum based on realistic explosion models, it 
is not certain that the models really give a good 
representation of material at velocities > 10, 000 
km s^^ since they emit little at the late-phases. 
Our explosion models contain material up to ^ 
15, 000 - 25, 000 km s'^ for fv = l- 1-6, above 
which the density drops very rapidly as a function 
of radius and therefore was not traced by the two- 
dimensional explosion calculations. However, the 



early-phase spectra suggest there is material up to 
~ 40, 000 km s~^ which carries a substantial frac- 
tion of the total kinetic energy (Nakamura et al. 
2001; see also Mazzah ct al. 2000). Although the 
origin of such material is not clear, its existence 
suggests that the total kinetic energy should be 
larger than that of the present two-dimensional 
model. 

In any case, the lower limit of the kinetic energy 
E51 > 10 is much larger than the typical value 
£51 ~ 1. Less energetic models do not fit the 
observed nebular spectra of SN 1998bw, especially 
at [Fell] 52OOA. 

Also, our models suggest that the ejecta is likely 
very clumpy and that nucleosynthetically differ- 
ent layers are mixed. This is necessary to explain 
the spectral evolution, especially the slow fading 
of [01] 6300A and Mgl] 4570A between +125'^ and 
-1-200'^. The amount of mixing is small and may 
well be explained in the context of Rayliegh- Taylor 
or shear instabilities between the 0-rich and the 
Fe-rich region (for Fe and Ca) and between the 
0-rich region and the C layer above it (for C). 
Such mixing processes are believed to take place 
in core-collapse supernovae (e.g., the earlier than 
expected detection of X-rays and 7-rays from SN 
1987A: Dotani et al. 1987; Sunyaev et al. 1987; 
Matz et al. 1988: see e.g., Kifonidis et al. 2000 for 
recent simulations). Also, aspherical explosions 
are suggested to boost the efficiency of the mixing 
(Nagataki, Shimizu, & Sato 1998). While we have 
introduced a rather large C abundance (10%) in 
the 0-rich region, both C and O arc hydrostatic 
burning products, and therefore additional mixing 
could also operate in hydrostatic evolution stages. 
This may be naturally explained by rotationally 
induced mixing (e.g., Iwamoto et al. 2005), since 
the large asphericity at the explosion implies that 
the progenitor star was a very rapid rotator. 

4.3. Implications for the Light Curve 

We derived a mass of ^^Ni, M{^'^Wi) ~ 0.4Mq 
in the two dimensional model. This is consistent 
with the results of previous work (e.g., Nakamura 
etal. 2001). The distribution of ^^Ni in our model 
is characterized by a large amount of ^^Ni in the 
high velocity outer region, and a small fraction in 
the inner high density region. The distribution 
will affect the shape of the light curve. Figure 10 
shows the evolution of the luminosity deposited 
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by 7-rays and positrons, i.e., optical light curve 
applicable after ~ 50 — 100 days. At least qualita- 
tively, the 2D models are favored over the original 
ID CO138E30 model, thanks to the presence of 
^^Ni at high velocities and the high density core 
at low velocities (Maeda et al. 2003a). Computa- 
tion of 2D light curves from the explosion to the 
late-phases on the basis of the two dimensional 
models will be another check of the validity of the 
models (see also Hoflich, Wheeler, & Wang 1999). 

4.4. Notes for Future Observations 

Our suggestion that hypernovae are aspherical 
explosions may further be able to be confirmed by 
future observations. First, we suggest to take late- 
phase near-infrared (NIR) spectra of supernovae 
with SN 1998bw-likc early-phase spectra, as has 
been done for SNe la (e.g., Hoflich et al. 2004). In 
the NIR, a [Fell] line is more isolated than in the 
optical band, making the effect of the orientation 
easier to see. 

It is also interesting to examine the [01] 6300A 
profile in other hypernovae and SNe Ib/c to inves- 
tigate the effect of the viewing angle. Recently, 
Kawabata et al. (2004) reported the detection of 
the double-peaked [01] 6300A in a spectrum of SN 
2003jd at ~ 1 year after the discovery. Mazzali et 
al. (2005) interpreted it as a SN 1998bw-like event, 
as suggested from early-phase spectroscopy, but 
viewed from the equatorial direction. Although 
the sample is small in number for hypernovae, late- 
phase spectra such as those published for SNe Ib/c 
by Matheson et al. (2001) would allow statistical 
studies to constrain the explosion energy and the 
asphericity of SNe Ib/c. 

4.5. Remaining Problems 

4.5.1. The Spectrum at +125^ 

One of the major advances in the present study 
is the time-dependent computations of the nebu- 
lar spectra of SN 1998bw. The model derived at 
the spectrum at -1-337'^ gives nice fits to the se- 
quence of the spectra after +200'^. However, at 
the first epoch -|-125'^, some deviation from the 
observations still exists especially in the luminos- 
ity of the [01] 6300, 6363A doublet, even though 
we have tried to fix this by introducing mixing and 
clumping. Others are the strong [Fell] and [Felll] 
between 4000A and 5000A and the strong Mgl] 



4570A and [01] 5577A in the model. They may 
simply be due to inappropriate treatment of the 
first epoch -1-125'* in our computation since at this 
epoch the ejecta density is still high and the neb- 
ular representation may not totally be a good ap- 
proximation. Temporal evolution of the intensity 
of the [01] 6300 A changed around +125^* (Patat et 
al. 2001), implying that at this epoch the ejecta 
are still in the transition from the photospheric 
phase to the nebular phase. For example, the ob- 
servationally strong 01 7800A implies that allowed 
transitions may be dominant in some wavelength 
ranges. 

Because of the possible presence of emission 
lines of some strong allowed Fe II transitions, the 
identification of the 5200A feature as the blend of 
forbidden lines only may be uncertain. Axelrod 
(1980) computed optical depths of the strongest 
allowed Fe II lines near 52OOA for a model SN 
la nebula at days 87 and 264 to be of the or- 
der of 10 and 0.1, respectively. Given the larger 
ratio Ek/M{'^'''Ki) in our models for SN 1998bw 
than for SNe la, the optical depths will be even 
smaller. Rough estimates suggest that a few lines 
may have optical depth of order unity at the first 
epoch (-1-125''). Therefore, at this epoch the con- 
tribution of the allowed Fe II lines and related ra- 
diation transfer effects, e.g., line scattering, may 
affect (but probably not dominate) the shape of 
the 52OOA feature. At the same time, the same 
estimate shows that the optical depths of the al- 
lowed Fe lines will be very small (of the order of 
0.1 or less) at -1-200'* and thereafter, justifying the 
assumption that the contribution of the allowed 
lines and radiation transfer effects are negligible. 
See also Maeda et al. (2002) for a detailed discus- 
sion of the identification of the 5200A feature as 
the forbidden lines. 

As for the Mgl] 4570A, the line cmissivity is 
quite sensitive to the treatment of the photoion- 
ization radiation field because of its low ioniza- 
tion threshold (Houck & Fransson 1996; Kozma 
& Fransson 1998a), while the ionization by UV 
photons is not included in our present spectrum 
synthesis calculations (See also, e.g., SoUerman et 
al. 2004 for the effect of photoionization). This 
could partly be a reason of our failure in fitting 
the temporal evolution of the Mgl] 4570A lumi- 
nosity, since the effect is expected to be stronger 
at earlier epochs. We plan to extend our code 
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to allow better treatment at such relatively early 
epochs. 

4-5.2. Peculiar [Call] Emission 

The final question we should answer in the 
future concerns the origin of the peculiar [Call] 
7300A emission. As mentioned in Appendix, the 
model spectra are too red. It is interesting that 
another hypernova, SN 2002ap, shows the [Call] 
7300A line exactly at the correct position (Kawa- 
bata et al. 2002; Leonard et al. 2002; Wang et al. 
2003; Foley et al. 2003). The blueshift of the line 
in SN 1998bw seems to be unique. 

The feature is a complex blend and therefore 
it is difficult to clariiy the reason of the failure 
in the model. With this caveat taken in mind, 
we speculate possibilities which may explain the 
discrepancy. The problem is possibly related to 
the distribution of Ca. There are at least two 
possibilities to reconcile the problem. First, the 
distribution of Ca may be similar to O, rather 
than Fe. It will then be challenging to the the- 
ory of explosive nucleosynthesis. Another possible 
interpretation would be that SN1998bw may have 
asymmetry even between the two hemispheres. If 
we look at the event from the hemisphere with the 
larger kinetic energy, the center of the Ca distribu- 
tion would be blue shifted. This should however 
also apply to [Fell] 5200A, which does not seem 
to show the shift. The significantly blended na- 
ture of [Fell] 5200A may wash that signal away. 
In this context, SN2002ap would be an explosion 
in which the degree of this asymmetry is small or 
the viewing orientation is relatively large. This is 
an interesting possibility, and we will pursue this 
issue in the future. The above two possibilities 
could be resolved by examining a number of hy- 
pernova late-phase spectra. In our interpretation, 
the sample contains hypernovae with various view- 
ing angles, then the two different Ca distributions 
- oxygen-like or iron-like - should give different 
statistics. 
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grant-in-Aid for Scientific Research and the 21st 
Century COE program of the Ministry of Educa- 
tion, Culture, Sports, Science and Technology in 
Japan. KM is Research Fellow of the Japan Soci- 
ety for the Promotion of Science (JSPS). 
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A. Emission Line Profiles 
A.l. [Fell] 5200 A 

Figures 11, 12, and 13 show how the shape of the [Fell] blend at 5200A depends on the degree of asphericity 
and the observer direction. The dependence on the kinetic energy can also be seen by comparing Figures 11 
(/y=1.6), 12 (/v=1.15), and 13 (/y=0.7). 

The dependencies can be summarized as follows. (1) For a given asphericity, a smaller viewing angle 6 
leads to a broader feature (with the obvious exception of the spherically symmetric models). (2) Larger 
asphericity leads to a broader [narrower] feature for small [large] 9. (3) Larger kinetic energy leads to a 
broader feature irrespective of the degree of asphericity and orientation. 

The above dependencies can be understood from the distribution of '"'^Ni (Fe) in Figure 5. The distribution 
is elongated in the polar direction in the aspherical models. Larger asphericity leads to higher velocities in 
the polar direction and smaller ones in the equatorial direction. If the energy is larger, the velocity in all 
directions becomes larger. These characters of the ^^Ni distribution explain the dependence of the [Fell] 
feature on various parameters. 

The observed broad 5200A emission can be explained if the asphericity is large, the energy is large, and 
the viewing angle is small {9 < 30°). For example, Figure 11 shows that for these very energetic models 
(/y=1.6), the synthetic [Fell] 5200A is as broad as the observed one for BP > 2 given 9 < 30°. For less 
energetic models, the criterion is tighter: For /y=1.15, BP > 8 is necessary to produce the broad feature. 
For fv=0-7, the synthetic [Fell] 5200A is never as broad as the observed one, even if the asphericity BP 
is very large. Note also that if both the energy and the asphericity are extremely large (e.g., /y=1.6 and 
BP=16), the synthetic feature is too broad. 

A.2. [OI] 6300, 6363A 

The dependence of the [01] 6300, 6363A doublet profile on various parameters is shown in Figures 14 - 
16. The observed sharply-peaked profile with extended wings again favors highly aspherical models (BP16 or 
BPS) viewed from a direction close to the pole {9 < 30°). In addition, the narrow [01] 6300A sets constraint 
on the model in a different way from the way the [Fell] 5200 A line does. For the less energetic models 
(/v=0.7), all the aspherical models (BP = 2-16) give nicely peaked profiles if viewed near the pole. For 
/y=1.15, larger BP (BP > 4) is required. Furthermore, for /y=1.6, the synthetic [01] 6300A is broader 
than the observations, even for small viewing angles. 

The sharply-peaked [01] profile is produced by the edge-on view of a disk-like oxygen distribution. The 
model predicts a double-peaked profile if viewed near the equator as shown in Figures 14 - 16 for large 9. A 
spherical explosion model like CO138E30 produces a flat-topped profile because of the central hole in oxygen 
distribution (at ~ 10, 000 km s~^ for the kinetic energy E^i = 30; Figs. 4a and 5), which does not fit well the 
observations, irrespective of the total kinetic energy. This is especially evident in the larger kinetic energy 
model {fv > 1) which produces very broad [01] 6300A emission. 

A.3. Mgl] 4570A 

The observed profile of Mgl] 4570A is sharply peaked, similar to [01] 6300A (Patat et al. 2001). The 
computation of nucleosynthesis in a supernova explosion gives a distribution of Mg similar to that of O. 
Therefore, the line profile shows characteristics similar to those of O. Figures 17 - 19 show the synthetic line 
profiles of Mgl] 4570A for various models and viewing angles. The fit to this line is less certain than that of 
[01] 6300A because of significant contributions of the [Fell] and [Felll] lines. To examine the characteristics 
of pure Mgl] 4570A, the contribution of this line is also shown in Figures 17 - 19. The shape of the line 
favors highly-aspherical models viewed at small 9, supporting the results derived from the [01] 6300A profile. 
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A.4. [Call] 7291, 7324A 

Figures 20 22 show synthetic spectra of the present models around 7300A compared to the observed 
one at 337 days after the B maximum. The feature is a blend of [Call] 729lA, 7324A, [Fell] 7155A, 7172A, 
7388A, and 7452A, [Coll] 754lA, and [Oil] 7322A. In a typical nebular condition, the [Call] lines (whose 
contributions are individually shown in Figures 20 - 22) probably dominate at the line center, while [Fell] 
(and [Coll]) fills up the wings. The [Oil] 7322A line is weak in the model, because most oxygen is neutral, 
as suggested by the strong [01] 6300A line. 

As shown in the figures, the models are redder than the observations. As mentioned above, the feature is 
a complex blend making the reason of the failure difficult to identify. For possible origins of the blueshift, 
see §4. 

The [Call] 7292A and 7324A in the models BP16 with small viewing angle (61 ~ 0°) show double peaked 

profiles because the distribution of Ca, as a product of explosive nucleosynthesis, follows closely that of Fe 
rather than that of O, which is a product o f hydrostatic burning and a fuel at the explosion (Fig. 5). 
A higher degree of asphericity (i.e., larger BP) leads to more aspherically distributed Ca. Therefore, the 
double-peaked profile of Ca for small viewing angle {6 0°) can be interpreted as two-blobs of Ca moving in 
opposite directions observed on the symmetry axis. We do not see such a profile in SN 1998bw, suggesting 
a slightly off-axis viewing angle. 

Because of the significant blend of the feature at 7300A, constraining the energy in the ejecta by fitting 
this feature is very uncertain. As long as only line width is concerned, less energetic models are preferred if 
6 is smaller (e.g., 6 = 0°), and more energetic ones if 6 is larger (e.g., 6 = 30°), since narrowing [broadening] 
the line by smaller [larger] energy must be compensated by broadening [narrowing] the line by smaller [larger] 

e. 
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Table 1 
Models 



Model'' Asphericity (BP) fv^ Ek/IO^^ ergs M(^«'Ni)/M0 Mq/Mq 



BP16 16 1.9 

16 1.75 

16 1.6 

16 1.45 

16 1.3 

16 1.15 

16 1.0 

16 0.85 

16 0.7 

BPS 8 1.9 

8 1.75 

8 1.6 

8 1.45 

8 1.3 

8 1.15 

8 1.0 

8 0.85 

8 0.7 

BP4 4 1.9 

4 1.75 

4 1.6 

4 1.45 

4 1.3 

4 1.15 

4 1.0 

4 0.85 

4 0.7 

BP2 2 1.9 

2 1.75 

2 1.6 

2 1.45 

2 1.3 

2 1.15 

2 1.0 

2 0.85 

2 0.7 

BPl 1 1.9 

1 1.75 

1 1.6 

1 1.45 

1 1.3 

1 1.15 

1 1.0 

1 0.85 

1 0.7 



33.1 


0, 


,79 


7.6 


26.5 


0, 


.67 


7.0 


20.7 


0, 


.63 


6.4 


15.9 


0, 


.60 


5.8 


11.5 


0, 


.52 


5.0 


8.2 


0, 


.45 


4.4 


5.9 


0, 


.34 


4.1 


3.8 


0, 


.27 


3.3 


2.3 


0, 


.20 


2.8 


34.5 


0, 


,75 


8.7 


27.5 


0, 


,68 


8.1 


21.3 


0, 


.63 


7.3 


15.3 


0, 


,55 


6.1 


11.5 


0, 


,48 


5.6 


8.4 





,37 


5.1 


5.9 


0, 


,32 


4.6 


4.0 


0, 


,25 


4.1 


2.0 


0, 


,22 


2.6 


32.5 


0, 


,63 


7.9 


25.9 


0, 


,58 


7.3 


20.2 


0, 


,53 


6.7 


15.3 


0, 


,45 


6.0 


11.2 


0, 


,42 


5.3 


7.9 





,36 


4.5 


5.7 


0, 


,30 


4.2 


3.7 


0, 


,25 


3.5 


2.1 


0, 


,19 


2.6 


36.9 


0, 


,55 


9.7 


28.9 


0, 


,51 


8.7 


22.1 


0, 


,47 


7.8 


16.3 


0, 


,43 


6.7 


12.2 


0, 


,36 


6.1 


8.8 


0, 


,33 


5.5 


6.2 


0, 


,27 


4.9 


4.1 


0, 


,24 


4.3 


2.9 


0, 


,22 


3.4 


37.5 


0, 


,54 


9.5 


29.7 


0, 


,49 


8.6 


23.1 


0, 


,47 


7.7 


16.7 


0, 


,43 


6.4 


12.7 





,41 


5.8 


9.5 


0, 


,39 


5.4 


6.6 


0, 


,33 


4.7 


4.1 


0, 


,30 


3.5 


2.5 


0, 


,22 


2.9 



*The models corresponds to the ones in Maeda et al (2002) as follows: BP16 (A), BPS 
(C), BP2 (E), BPl (F: a spherical model). The models BP4 are newly calculated in the 

way same with Maeda et al. (2002). 

Velocities at any points in the ejecta are multiplied by a factor of fv- fv — 1 for 
the original models in Maeda et al. (2002). 
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Fig. 1. — Synthetic spectra (black) for the original CO138E30 model as compared with the spectra of SN 
1998bw (gray). For presentation, both the synthetic and observed spectra are multiplied by an arbitrary 
constant (consti) and then vertically shifted upward by another factor (const2). The amounts (consti, const2 
in units of lO-^^ gj.g g-i ^^-2 ^-1^ follows: +125'^ (0.1, 1.8), +200^^ (0.1, 1.2), +337'^ (0.5, 0.6), and 

+376<i (1.0, 0.0). The mass of ^^^Ni is O.8M0. 
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Fig. 2. — Synthetic spectra for the modified CO138E30 model. See the caption of Fig. 1. The mass of ^^Ni 
is O.61M0. 
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Fig. 3. Density distribution at 1 day after the explosion. The original CO138E30 (gray) and the modified 
model (black) are shown. High velocity materials at > 30, 000 km s~^ are not included in the present study, 
since they contribute virtually nothing on late phase emission. 
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Fig. 4. — Mass fractions of some selected elements and ^^Ni in the velocity space for (a) the original 
CO138E30 and for (b) the modified model. 
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Fig. 5.— Structure of the ejecta for Models BP16, BPS, BP4, BP2, and BPl (from left to right) in the 
ejecta velocity space at a reference time (at 100 seconds). The parameter fy = 1. The axes are r and z 
components of velocities scaled at 15, 000 km s~^. For each model, shown here are density, mass fraction of 
^^Ni, Ca, and O (from top to bottom). 
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Fig. 6. - Examples of synthetic spectra (black lines) as compared with the spectrum of SN 1998bw at 337 
days after the maximum light (gray lines). Shown here are Models BPS, BP2, and BPl with (a) /v=1.15 
and e = 30° and with (b) /y=1.6 and 6 = 30°. The spectra of BPS and BP2 are shifted vertically for 
presentation (+2 and +1 xlO~^^ erg cm"^ A~'^) for BPS and BP2, respectively). 
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Fig. 7. — Acceptable models for 1998bw. (a) 6 = 0° and (b) 6 = 30°. The figure is obtained by comparing 
model spectra with the observed one of SN 1998bw at 337 days after the B maximum. The regions surrounded 
by solid, dashed, and dotted lines contain models which give an acceptable fit to the line profiles around 
5200 A ([Fell]), 6300 A ([01]), and 7300A ([Call]), respectively. See §3.4 for details. 
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Fig. 8. — Temporal evolution of the synthetic spectra for 2D models as compared with the spectra of SN 
1998bw. Shown here are Model BPS with /v'=1.15 (a) and /y=1.6 (b). The orientation is 6* = 30°. For 
presentation, both the synthetic and the observed spectra are multiplied by an arbitrary constant (consti) 
and then vertically shifted upward by another factor (const2). The amounts (consti, const2 in units of 10~^^ 
erg s-i cm-2 A"!) are as follows: +125^ (0.1, 1.8), +200'^ (0.1, 1.2), +337'^ (0.5, 0.6), and +376^^ (1.0, 0.0). 
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Fig. 9. Sequence of the spectra of the "modified" two-dimensional models (black solid lines). Shown here 
are Model BPS with /v=1.15 (a) and /y=1.6 (b). The orientation is 6 = 30°. The spectra of SN 1998bw 
(gray lines) and the modified "one dimensional" CO138E30 model (black dotted line) are also shown. See 
also the caption of Figure 8. The mass of ""'^'Ni in the modified models is O.4OM0 for /y=1.15 (a) and O.56M0 
for /y=1.6 (b), respectively. 
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Fig. 10. — The bolometric light curve of SN 1998bw (gray points; Patat et al. 2001) arc compared with the 
synthetic curve from the ID original CO138E30 (thin lines) and that from the 2D BPS model with /v'=1.15 
(thick solid). The mass of ^^Ni in the 2D model is G.STM© (Table 1). In the ID model, the mass of ^''Ni 
is varied to iUustratc the difference with the 2D model. M(^^Ni) is O.64M0 (thin sohd) and O.4IM0 (thin 
dotted) for the ID model. 
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Fig. 11. — The synthetic spectra of the 2D models around 5200 A (black lines) as compared with the observed 
one of SN 1998bw at 337 days after the B maximum (gray lines). The synthetic spectra are mainly contributed 
by [Fell]. The models are those with /y = 1.6. Shown here are BP16, 8, 4, 2, 1 (from left to right) and 
6* = 0, 30, 60, 90° (from top to bottom). 
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Fig. 12. — The same with Fig. 11 ([Fell] 5200A), except the model parameter fv = 1.15. 



28 



5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 




5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 




5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 

2 I ' ' ' ' ' ' I 2 I ' I ' ' T ' ' ' ■ ' ■ I 2 I • I I 2 ' 




5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 5000 5200 5400 



Fig. 13. — The same with Fig. 11 ([Fell] 5200A), except the model parameter fv = 0.7. 
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Fig. 14. — The synthetic spectra of the 2D models around 6300A (black lines) as compared with the observed 
one of SN 1998bw at 337 days after the B maximum (gray lines). The synthetic spectra are predominantly 
contributed by [01] 6300A and 6363A. The models are those with fv = 1.6. See the caption of Fig. 11 for 
the meanings of each panel. 
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Fig. 15.— The same with Fig. 14 ([01] 6300A and 6363A), except the model parameter /y = 1.15. 
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Fig. 16.— The same with Fig. 14 ([01] 6300A and 6363A), except the model parameter fv = 0.7. 
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Fig. 17. — The synthetic spectra of the 2D model around 4600A (black thick lines) as compared with the 
observed one of SN 1998bw at 337 days after the B maximum (gray lines). The synthetic spectra are mainly 
contributed by Mgi] 4570A and forest of [Fell] and [Felll]. The models arc those with fy = 1-6. See the 
caption of Fig. 11 for the meanings of each panel. The contribution of Mgl] 4570A is also shown (thin lines). 



33 





4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 




4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 




4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 

2 I I I I I I I I I I I 2 I ' I ' ' ' I ' ' ' I ' I 2 I ' I ' ' ' I ' ' ' I ' I 2 I ' I ' ' ' I ' ' ' I ' 




4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 4400 4600 4800 



Fig. 18.— The same with Fig. 17 (Mgl] 4570A and [Fell], [Felll]), except the model parameter fy = 1.15. 
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Fig. 19.— The same with Fig. 17 (Mgl] 4570A and [Fell], [Felll]), except the model parameter fv = 0.7. 
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Fig. 20. — The synthetic spectra of the 2D models around 7300A (black thick lines) as compared with the 
observed one of SN 1998bw at 337 days after the B maximum (gray lines). The synthetic spectra are mainly 
contributed by [Call] 729lA, 7324A, and [Fell] 7155A, 7172A, 7388A, and 7452A. The models are those 
with fv = 1.6. See the caption of Fig. 11 for the meanings of each panel. The contributions of [Call] 729lA 
and 7324A are also shown (thin lines). 
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Fig. 21.— The same with Fig. 20 ([Call] 729lA, 7323A and [Fell]), except the model parameter fy = 1.15. 
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Fig. 22.— The same with Fig. 20 ([Call] 729lA, 7323A and [Fell]), except the model parameter fy = 0.7. 
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